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Abstract 

It is shown that a linear separation relations are fundamental objects for integration by quadra- 
tures of Stackel separable Liouville integrable systems (the so-called Stackel systems) . These relations 
are further employed for the classification of Stackel systems. Moreover, we prove that any Stackel 
separable Liouville integrable system can be lifted to a bi-Hamiltonian system of Gel'fand-Zakharevich 
type. In conjunction with other known result this implies that the existence of bi-Hamiltonian rep- 
resentation of Liouville integrable systems is a necessary condition for Stackel separability. 

1 Introduction 

The Hamilton- Jacobi (HJ) theory seems to be one of the most powerful methods of integration by 
quadratures for a wide class of systems described by nonlinear ordinary differential equations, with a 
long history as a part of analytical mechanics. The theory in question is closely related to the Liouville 
integrable Hamiltonian systems. The milestones of this theory include the works of Stackel, Levi-Civita, 
Eisenhart, Woodhouse, Kalnins, Miller, Benenti and others. The majority of results was obtained for a 
very special class of integrable systems, important from the physical point of view, namely for the systems 
with quadratic in momenta first integrals. 

The first efficient construction of the separation variables for dynamical systems was discovered by 
Sklyanin [1]. He adapted the methods of soliton theory, i.e., the Lax representation and r-matrix theory 
for systematic derivation of separation coordinates. In this approach the integrals of motion in involution 
appear as coefficients of characteristic equation (spectral curve) of the Lax matrix. This method was 
successfully applied for separating variables in many integrable systems [l]-[8]. 

Recently, a modern geometric theory of separability on bi-Poisson manifolds was developed [9]-[15]. 
This theory is closely related to the so-called Gel'fand-Zakharevich (GZ) bi-Hamiltonian systems [16], [17]. 
The theory in question includes Liouville integrable systems with integrals of motion being functions 
quadratic in momenta as a very special case. In this approach the constants of motion are closely related 
to the so-called separation curve which is intimately related to the Stackel separation relations. The 
separation curve arising in the geometric approach is closely related to its counterpart in the r-matrix 
approach. In fact, these curves are identical for linear r-matrix and related by exponentiation of momenta 
in the spectral curve for dynamical (quadratic) r-matrix [6], [14]. 

In the present paper we develop in a systematic fashion a separability theory of the Liouville integrable 
systems which are of the GZ type, including as a special case the class of systems with quadratic in 
momenta first integrals. First of all, we treat Stackel separable systems according to the form of separation 
relations and make some observations related to their classification. Then we construct a quasi-bi- 
Hamiltonian representation of Stackel systems on 2n dimensional phase space and lift them to the related 
GZ bi-Hamiltonian systems on the extended 2n + k dimensional phase space. This result proves that 
bi-Hamiltonian property is common to all classes of the Stackel systems considered. In other words, we 
prove that the existence of bi-Hamiltonian representation for Liouville integrable systems is a necessary 
condition for their Stackel separability, i.e. for these systems for which separation relations are linear 
in all constants of motion which are in involution. Finally let us mention that up to now such a proof 



was available only for a distinguished class of the so-called Benenti systems [22], where k = 1 and all 
constants of motion are quadratic in momenta. 



2 Separable Stackel systems 



Consider a Liouville integrable system on a 2n-dimcnsional phase space M. Thus, we have M 3 u = 
(qi, . . . ,q n ,Pi, ■ ■ ■ ,Pn) T an d there are n functions Hi(q,p) in involution with respect to the canonical 
Poisson tensor n 

{HijHjjn = ir(dHi,dHj) = (dH^irdH^) >= 0, i,j = 1, . . . ,n, 

where (•, •) is the standard pairing oiTM and T*M. Canonicity of 7r means that the only nonzero Poisson 
brackets among the coordinates are {qi,Pj}v — hj- The functions Hi generate n Hamiltonian dynamic 
systems 

u u = tt dH l = X Hi , i = l,...,n, (1) 

where Xji i are called the Hamiltonian vector fields. 

The Hamilton- Jacobi (HJ) method for solving (1) essentially amounts to the linearization of the latter 
via a canonical transformation 

(q,p) -> (b, a), a.^H,, i = l,...,n. (2) 

In order to find the conjugate coordinates bi it is necessary to construct a generating function W(q, a) of 
the transformation (2) such that 

b- = 9W _ dw 

The function W(q,a) is a complete integral of the associated Hamilton- Jacobi equations 

H i (q 1 ,...,q n ,^,...,^-)=a i , t = l,...,n. (3) 
In the (6, o) representation the ^-dynamics is trivial: 

( a j)u = 0, (bj)u = Sij, 

whence 

6 j(9»°) = ^ : =*j+ c j' j = l,-..,n, (4) 

where Cj are arbitrary constants. 

Equations (4) provide implicit solutions for (1). Solving them for qj is known as the inverse Jacobi 
problem. The reconstruction in explicit form of trajectories qj = qj(U) is in itself a highly nontrivial 
problem from algebraic geometry which is beyond the scope of this paper. 

The main difficulty in applying the above method to a given Liouville integrable system in given 
canonical coordinates (q,p) consists in solving the system (3) for W. In general this is a hopeless task, 
as (3) is a system of nonlinear coupled partial differential equations. In essence, the only hitherto known 
way of overcoming this difficulty is to find distinguished canonical coordinates, denoted here by (A, fi), 
for which there exist n relations 

tp i (\ i ,p i ;a 1 ,...,a n ) = 0, i = 1, . . . , n, a, E R, dct |^ ^0, (5) 

such that each of these relations involves only a single pair of canonical coordinates [1] . The determinant 
condition in (5) means that we can solve the equations (5) for a, and express Oj in the form otj = i?j(A, /i), 
i = 1, . . . ,n. 

If the functions Wi(Xi,a) are solutions of a system of n decoupled ODEs obtained from (5) by 
substituting ^ — dWi ( Xi ' a '> 



Vi (*i,IH= dW ^'" ) ,ai,...,o n ) ^0, i = l,...,n, (6) 
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then the function 

W(X,a) = V" W^a) 
z — *i=i 

is an additivcly separable solution of all the equations (6), and simultaneously it is a solution of all 
Hamilton- Jacobi equations (3) because solving (5) to the form ai = i?j(A, /i) is a purely algebraic op- 
eration. The Hamiltonian functions Hi Poisson commute since the constructed function W(X, a) is a 
generating function for the canonical transformation (A, ii) — > (6, a). The distinguished coordinates (A, ii) 
for which the original Hamilton- Jacobi equations (3) are equivalent to a set of separation relations (6) 
are called the separation coordinates. 

Of course, the original Jacobi formulation of the method was a bit different from the one presented 
above, and was made for a particular class of Hamiltonians, nevertheless it contained all important ideas 
of the method. Jacobi himself doubted whether there exists a systematic method for construction of 
separation coordinates. Indeed, for many decades of development of separability theory, the method did 
not exist. Only recently, at the end of the 20th century, after more then one hundred years of efforts, two 
different constructive methods were suggested, the first related to the Lax representation and the second 
related to the bi-Hamiltonian representation for a given integrable system. 

We would like to stress that all results of the present paper are derived directly from the separation 
relations (5), thus confirming their fundamental role in the modern separability theory. 

In what follows we restrict ourselves to considering a special case of (5) when all separation relations 
are affine in H^. 

n 

^S^{\i,^i)H k =^i{\i,Hi), i = l,...,n, (7) 
fc=i 

where S k and ipi are arbitrary smooth functions of their arguments. The relations (7) are called the 
generalized Stackel separation relations and the related dynamical systems are called the Stdckel separable 
ones. The matrix S = (S k ) will be called a generalized Stdckel matrix. The reason behind this name is 
the fact that the conditions (7) with S k being /x-independent and ipi being quadratic in momenta /i are 
equivalent to the original Stackel conditions for separability of Hamiltonians Hi. To recover the explicit 
Stackel form of the Hamiltonians it suffices to solve the linear system (7) with respect to Hi. 

Although the restriction of linearity appears to be very strong, for all known separable systems (at 
least to the knowledge of the author), the general separation conditions can be reduced to the form (7) 
upon suitable choice of integrals of motion Hi. The possible explanation of this fact is that we simply 
have no mathematical tools for effective construction of separation coordinates for non-Stackel separable 
systems, so that part of separability theory is yet terra incognita. 

Let us come back to the Stackel case. If in (7) we further have S k (Xi, /x,) = S k (Xi, /Xj) and ipi(\i, Hi) — 
ip(Xi,Hi) then the separation conditions can be represented by n copies of the curve 

n 

Y,S k {\,li)H k =i>(\n) (8) 

k=l 

in (A, pi) plane, called a separation curve. The copies in question are obtained by setting A = Aj and 
/i = fii for i = 1, . . . ,n. 

Remark. There is an important special case when (8) is an arbitrary nonsingular compact Ricmann 
surface T, i.e., when S k (X,[i) and ^>(A,/x) are polynomials of A and /x of certain specific form. Then one 
can find the genus of this curve and basic holomorphic differentials in a standard fashion and the Jacobi 
inversion problem (4) can be equivalcntly expressed by the Abel map of the Ricmann surface T into its 
Jacobi variety and solved in the language of Riemann theta functions (see [18] and references therein). 

From now on we will consider Stackel separable systems with separation relations of the most general 
form (7). For reasons to be explained in the next section, we collect the terms from the l.h.s. of (7) as 
follows: 

m 

^^(A l ,^)^ (fc) (A 4 )=^(A l , M4 ), i = l,...,n, (9) 
fe=i 

where 

H<- k \X)=J2 A " fc ' »i + • ■ ■ + n m = n 

i=l 
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and impose the normalization <^™(Aj,/ij) = 1. 

As the separation relations (7) play the fundamental role in the Hamilton- Jacobi theory, it is natural 
to employ them for classification of Stackel systems. The form of separation relations (9) allows us 
to classify the associated Stackel systems. Actually, any given class of Stackel separable systems can 
be represented by a fixed Stackel matrix S and the functions ip. The matrix £ is uniquely defined by 
to vectors ip k = (<pi, . . . , fn) T , k = 1, . . . , to, and the partition (m, . . . , n m ) of n. Note that in our 
normalization we have tp m = (1, . . . , 1) T . 

For example, the most intensively studied systems in the 20th century, those related to one-particle 
separable dynamics on Ricmannian manifolds with flat or constant curvature metrics, belong to the 
simplest class with m = 1 and the functions ^ being quadratic in the momenta 

n 

J2 H oK~ j = \fi^i)l4+!i^il t=l,...,n. (10) 
i=i 

This class, which will hereinafter be referred to as the Benenti class, includes systems generated by 
conformal Killing tensors [19]-[22], as well as bi-cofactor systems, generated by a pair of conformal Killing 
tensors [23]- [27]. Here the functions /j define the Stackel metric while the functions 7, define a separable 
potential. When fa = /(Aj) and / is a polynomial of order not higher then n + 1, then the associated 
Stackel metric is of constant curvature. 

Another class of separable systems also has to = 1 but the functions are now exponential in the 
momenta 

n 

Hj\™~ 3 = exp(aft) + exp(-fyuj) + 7,(Ai), i = l,...,n, 

3=1 

where 7, define a separable potential. This class includes such systems as the periodic Toda lattice [13], 
the KdV dressing chain [14], the Ruijsenaars-Schneider system [11] and others. 

We also know some particular examples from the classes with to > 1. For instance, stationary flows 
of the Boussinesq hierarchy belong to the class with to = 2, m = 1, ni = n — 1, f} = m [12], [11]. 
Dynamical system on loop algebra s[(3) belongs to the class with to = 2, m = 2, n 2 — 4, <p\ = fii [15]. 
In both cases the functions tpi are cubic in the momenta, so these separation relations belong to the 
following class: 

ni n 2 

Hi Hf ] \T~ 3 + E H f X ?~° = \f^i)t4 + mi(Ai) + 72(A 4 ), i = 1, . . . , n, (11) 

i=i i=i 

where fi'fi and 72 give rise to the separable potentials. 

Finally, systems from the classes with 1 < to < n, ip\ — x\ ak \ «fe £ N and with ipi quadratic in the 
momenta, i.e. 

m 

y £\^ ) H( k \\ l ) = lf t (\ i )tf+ji(\i), i = l,...,n, (12) 
k=i 

were constructed in [31] and [32]. 

3 Bi-Hamiltonian property of Stackel systems 

We start this section with a few definitions important for further considerations. As the Hamiltonian 
formalism is of tensorial type, there is no need to restrict ourselves to nondegenerate canonical repre- 
sentation of Hamiltonian vector fields. Given a manifold A4, a Poisson operator ir on M. is a bivector 
(second order contravariant tensor field) with vanishing Schouten bracket 

[ir,ir] s = 0. 

Then the bracket 

{h,hh ■■= (dhMh), fuh e c°°(M), 
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is the Lie bracket, i.e., it is skew-symmetric and satisfies the Jacobi identity. A function c : M. — > M is 
called the Casimir function of the Poisson operator it if for an arbitrary function / : A4 — > K we have 
{/, c}^ = (or, equivalently, if 7rdc = 0). A linear combination tt\ = m — Xttq (A G R) of two Poisson 
operators ttq and 7Ti is called a Poisson pencil if the operator 7Ta is Poisson for any value of the parameter 
A, i.e., when [7r ,7ri]s = 0. In this case we say that 7r and 7Ti are compatible. Given a Poisson pencil 
n\ = m — Xn we can often construct a sequence of vector fields Xi on M that have two Hamiltonian 
representations (the so-called bi- Hamiltonian chain) 

Xi = TTxdhi = Tr dh i+1 , (13) 

where hi G C°°{M) are called the Hamiltonians of the chain (13) and where i is a discrete index. This 
sequence of vector fields may or may not terminate in zero depending on the existence of the Casimir 
functions for the pencil. 

Consider a bi-Poisson manifold (M, ttq, tti) of dim M = In + m, where 7To, 7Ti is a pair of compatible 
Poisson tensors of rank In. We further assume that the Poisson pencil tt\ admits m Casimir functions 
which are polynomial in the pencil parameter A and have the form 



h^(\)=J2\ n i- i h < ?\ j = l,...,m, (14) 

i=0 



(7) 

so that m + ■ ■ ■ + n m = n and h\ arc functionally independent. The collection of n bi-Hamiltonian 
vector fields 

nxdh^iX) = ^ Xl j) = mdh^ = wodh^, i=l,..., nj , j = l,...,m, (15) 

is called the Gel'fand-Zakharevich (GZ) system of the bi-Poisson manifold M. Notice that each chain 
starts from a Casimir of ttq and terminates with a Casimir of tt\. Moreover, all h\^ pairwise commute 
with respect to both Poisson structures 

xPihP) = (dh\ k \n dh^ = {dhf\^dh^) = {hflh^U = {h^,h^U = 0. 

In the following section we prove that an arbitrary Stackel system on the phase space M with separa- 
tion conditions given by (9) can be lifted to a GZ bi-Hamiltonian system on the extended phase space M . 

As recently proved in [15], the Stackel Hamiltonians from separation relations (7) admit the following 
quasi-bi-Hamiltonian representation 

n 

n 1 dH i =^F ij n dH j , * = l,...,n, (16) 



where n is a canonical Poisson tensor 

n n 



/„ 

-In 



I n is an n x n unit matrix, III is a noncanonical Poisson tensor of the form 

III = ( ^ V J , A„ = diag(Ai,...,A„), 



-A n 

compatible with n , and the control matrix F has the form 

F = (S~ 1 A n S), (17) 

where S is the associated Stackel matrix. 

To have a better insight into the functions Fij, we will find another representation for the entries 
Fij = (S~ 1 A n S)ij of F. To this end consider a system of n linear equations for Vk, k = 1, . . . , n; 

n n 

J2sh^»i)Vk = J2 X i S i( X "^ a l> i = l,...,n, (18) 
fc=i i=i 
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where a,i,i = 1, . . . , n are some parameters. The solution of this system has the form 

v - V* - det(S™) 

p=l 

where S^p) is the matrix S 1 with the r-th column replaced by (XiSf (Ai/Ui), . . . , A„S'P(A„/i„)) T , the string 
of coefficients at the parameter a p . On the other hand, as V = (V\, . . . , V n ) T and {a\, . . . , a n ) T , system 
(18) can be written in the matrix form as 

JV = A n Ja =^> V = J~ 1 A n J a = a a, 

where a>ij = (J~ 1 A n J) i j. Comparing this result with (17) and (19) we find 

Fa = (J-'AnJh = (20) 

Now, the important question is which entries F^ are nonzero when the separation relations take the 
form (7). In other words, we want to know for which i,j det(S^) ^ 0, i.e., the matrix has no 
linearly dependent columns. 

To answer this question, we first rewrite the quasi-bi-Hamiltonian chain (16) in the equivalent form 

m ni 

n x dH\ k) =J2J2 F ?jn dH«\ k = l,...,m, i = l,...,n k (21) 

1=1 j = l 

adapted to the separation relations written in the form (9). Then a simple inspection shows that 

Ftf+i = 1, F^= F k ' 1 ± 0. 
Hence, the quasi-bi-Hamiltonian representation (21) takes the form 



Yl.dH^ = H dffW + J2 Ho dH® , = 0, (22) 



1=1 

where 



k>l = detsj fc -° 
detS ' 

and 5} fe '^ is the Stackel matrix 5 with («! + •■• + "fc-i + i)-th column replaced by (^A"', . . . , ip l n A™ l ) T . 

In the rest of this section we show that the representation (22) can be lifted to a GZ bi-Hamiltonian 
form. First, we extend the 2n dimensional phase space M to M. = M x M m with additional coordinates 
Ci, i — 1, . . . , to, on M m . Then, we extend the Hamiltonians as follows: 



H\ k) (q, p) - ftj fc) (9, P, c) - tff > (9, p) - («> P) c ' • ( 23 ) 

From (18)-(20)we infer that the separation relations for read 

m 

^^(A 4 ,^)^ (fc) (A 4 ) =Vi(Ai,w), i = l,...,n, (24) 
fe=i 

where 

/ l W(A)=^A" fc " i / l f ) , 4 fe) =c fc , m + --- + n m =n. 

i=0 

Moreover, for the functions F k ' 1 we have the same quasi-bi-Hamiltonian representation as for H- k ^: 

m 

U,dF k ' 1 = n dF k +[ + F? ,r Ho dFl' 1 , (25) 



r=l 
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This was proved for arbitrary Stackel systems in [15]. 

Denote the push- forwards of the Poisson tensors n and III to M by ttq and tt\d- Both tto and ttid 
are degenerate and possess common Casimirs c i7 i — 1, . . . , to. We have 



n 












ni 












(26) 



Relations (22)- (26) imply that on M we have a quasi-bi-Hamiltonian representation with respect to the 
Poisson tensors ttq and -kid of the form 

m 

n 1D dh^ ^^dh^+^F^nodh^, F^ l = -6 kl , h™ +1 = 0. (27) 
i=i 

Now introduce the bivector 

m 

7T1 -TTID+^X^ AZ k , 
k=l 

where 

First we show that the bivector tti is Poisson. Using the properties of the Schouten bracket we have 



[tti, 7Ti] s = 2 ^ A L^ottid + 2 ^[X^, Z,-] A Zi A x[ 3 \ (28) 



where means the Lie derivative in the direction of X, and [•, •] is the commutator of vector fields. 
Now, let us prove that 

L x p Tl 1D =J2*o dF r / A X[ L) . (29) 
1 l 

From (27) we have 

I 

r 

From the Poisson property of nm it follows that 

= Y,(Fnt L x{^D - 7T1D dF% A X[ l) ) 

I 

= E V 10 - E F - ^ A 
= E V 1 ^ - E F - E ^ A 



= E ^ r ^ V )7riD - E *° dF ^ 1 A . 



hence (29) is satisfied. On the other hand, 

[X®,Z j ]=n dFi*, 

so (28) becomes 

[tti , m] s = 2 ^ Zi A tto dF[' J A x[ j) + 2 ^ tt dF^ A Z t A x[ j) = 0, 



i,3 »>J 
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and thus 7Ti is Poisson. 

Moreover, the Poisson bivectors ir and m are compatible as 

(i) 



[7r ,7ri] 5 = ^(Zi A L X ^U - X[ l) A L Zi n ) 



0. 



(k) 

Finally, the vector fields X!> form bi-Hamiltonian chains with respect to 7ro,7ri. Indeed, we have 



1=1 



as 



Quite obviously, 



a zi)dh\ k) = -j2 F i' 1 x i l) = - E n " dfe i°- 
i=i i=i i=i 



m dh%> = o, x[ l) = in dci = ni d/l^ , 

so h^(X) are polynomial in A Casimir functions of the Poisson pencil it\ = m — Xttq- 



4 Examples 

Here we illustrate the above results with three examples of separable systems with three degrees of 
freedom. Two of them are classical Stackel systems with separation relations quadratic in the momenta 
while the third example has separation relations cubic in the momenta. 
Example 1. 

Consider the separation relations on a six-dimensional phase space given by the following bare (potential- 
free) separation curve 

ffiA 2 + JT 2 A + JT 3 = |m 2 

from the class (10). This curve corresponds to geodesic motion for a classical Stackel system (of Benenti 
type). The transformation (A, fi) — ► (q,p) to the flat coordinates of associated metric follows from the 
point transformation 

ai(q) = qi = X 1 + X 2 + A 3 , a 2 (q) = \q\ + \q 2 = A 1 A 2 + A 1 A 3 + A 2 A 3 , a 3 (q) = \q x q 2 + \q 3 = A 1 A 2 A 3 . 
In the flat coordinates the Hamiltonians take the form 
Hi =pip 3 + \p\, 

h 2 = \qzpl - \ qipI + \q2P2Pz - \pipi - \qiPlVi, 

H 3 = \q\v\ + \q\p\ + \v\ + \qiviP2 + \q2p1pz - \q\qipipz - \qzviP-s 



and admit a quasi-bi-Hamiltonian representation (16) with the operators n and IT of the form 



n 



o / 3 
-/ 3 o 



n 1 = i 



/ 











qi 


-1 





\ 













92 





-1 















293 


92 


9i 






-5i 


-92 







P2 


P3 






1 





-92 


~P2 










V 





1 


-9i 


~P3 








/ 



(30) 



(31) 
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and the control matrix 



qi 1 

.12 _ 12 
4*1 2' 

ki«2 + 3«s o o 



\qi o i 

On the extended phase space of dimension seven, with an additional coordinate c, we have the extended 
Hamiltonians (23) 



h = c, 



hi =PiP3 + \v\ -cqi, 



>q3P 3 



h 3 = \q\v\ + \q\v\ + gPi + \q1P1P2 + \q2V\Vz - \q\q2V2Vi - \q-iV2Vz - {\q\q2 + \qs)c 

They form a bi-Hamiltonian chain 



ir dh 


= 




TTodhl 


= x 1 


= iridho 




= x 2 


— -K\dhi 


ir dh 3 


= X 3 


= -Kxdh-i 







= TT\dh 3 ^ 



with the Poisson operators wq and tt\, where 



7T = 



n 












7Ti = 





Xy ' 








The separation curve for the extended system takes the form 

cA 3 + /11A 2 + h 2 X + h 3 = 



Example 2. 

Consider now separation relations on a six-dimensional phase space given by the following bare separation 
curve 

i?iA 3 + H 2 X 2 + H 3 = 
from the class (12). When written using the notation (9), this curve takes the form 



X^H^X + H^+H^ 



and again represents geodesic motion for a classical Stackel system (this time of non-Benenti type) . Using 
the coordinates, the Hamiltonians, and the functions Ui from the previous example we find that 



#1 


= 


= -±H 2 , 

172 ' 


H 2 


= hP 


= ^1-^2 


H 3 




= H 3 -%H 2 



and thus we see that the Hamiltonians Hi are related to Hi through the so-called generalized Stackel 
transform (see [32] for further details on the latter). One can show that the metric tensor associated to 
the Hamiltonian Hi is not flat anymore. 

The Hamiltonians Hi form a quasi-bi-Hamiltonian chain (16) with the Poisson tensors (30), (31) and 
the control matrix 

CTi-^i 1 - — 

(72 (72 

-a 2 + 

2 <T2 
a 3 Q £3 

(72 72 



9 



On the extended phase space of dimension eight, with additional coordinates Ci,c 2 , we have the 
extended Hamiltonians (23) 

= ci, 
tf> = c 2 , 



They form two bi-Hamiltonian chains 
nodho = 



ir dh^ 
irodh^ 



X 
X 



(i) 
l 

(i) 



,(1) 



Tr dh y 



(2) 







TT dh[ 2) = X[ 



Tridh\ 
= Trid/i^ 

with the Poisson operators ttq and 7Ti of the form 



(2) 



iTidh^P 



(2) 





/ n 





TTO = j 










v o 



7T1 



ni 



-(*! 1} ) T 
V -(*i (2) ) T 



= TTld/li' 



The separation curve for the extended system takes the form 

A 2 ( Cl A 2 + h^X + + c 2 A + /if } = |m 2 - 

Example 3. 

Consider separation relations on a six-dimensional phase space given by the following bare separation 
curve cubic in the momenta 

,3 



from the class (11). The transformation (A, fi) — ► (<Z,p) to new canonical coordinates in which all Hamil- 
tonians are of a polynomial form is obtained from the following two transformations: 



ui 
u 2 

Vl 

v 2 



3q 2 - 3q 3 , 

-qip 2 - q\Pz + 3g 2 + 5^ - 69293, 

-ql - 9qfq 3 + qiq 3 p 2 + qiq 3 p 3 - 2JQ1Q2 + q\p\ + 3q- 2 gf, 

1 

Qi' 
3<? 2 - 2q 3 

1 

Qi 



2 93 
V 3 = p 3 + aP2 

6 Ql 



c,Q2q 3 . 2 

3 Aqf, 

Qi 



and 



ui = A1+A2 + A3, 

u 2 = A1A2 + A1A3 + A2A3, 
A1A2A3, 



"3 
Mi 



A, + -y 2 Ai + t) 3 , 



1,2,3. 
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In the ((7,p)-coordinates the Hamiltonians take the form 

# = P2P3 + \vl + P% - 7q 2 iP3 - ^q\p2 - 3<Z2Pi + 18gi?2 + 13«i + 12535i52, 
h[ 2) = I2q\q 2 + 8qfq 3 - 2q^p 1 + (-6qiq 2 - 4<2i<7 3 )jj 3 + pip 3 , 

tr( 2 ) 1 2,12 23 22422 2 10 2 

H 2 = + 3P2P3 + 4?P2 ~ QiP 3 ~ §5iP 2 - 92P1P2 - 5lPl - fqiP3P2 

+(53 - 3<7 2 )piP3 + (21<??<72 + 6<? 3 5i>i + (4<?3<Zi<?2 + 6gi<^ + f q\)p 2 
+ (7q$ + I8q iq i + 6g 3 5l52 - ^q^ps - 8qfq 2 - 72q 3 qfq 2 - 90qfq 2 2 - I2qf. 

They form a quasi-bi-Hamiltonian chain (16) with the non-canonical Poisson operator 



IIi = 



/ 











-53 


3?i 


252 


\ 










-\qi 


A 


3^2 - 53 


-52 









g5i 





2q\ 





-53 






-53 


-A 


-2q\ 





B 


C 






l-3gi 


-3<72 + 53 





—B 





-24q 2 




V 


2gi 


52 


93 


-C 


245? 





/ 



where A = —\p2 + \pz — 3(7?, B = 54qiq 2 + 24qiq 3 — 3pi, C = —2Aqiq 2 — 12qiq 3 + p\ and the control 
matrix 

(-53 -5i \ 

-\p2 + ql -2<Z3 + 3g 2 1 • 

5535i + 6<7?(72 - 5iPi - 3-53P2 -45i - 53 + 3 5253 + §5iP2 + 5iP3 / 

On the extended phase space of dimension eight, with additional coordinates ci,C2, the extended 
Hamiltonians (23) are 

4 1} = c ii 

= + q 3Cl + qi c 2 , 

1,(2) 

K = c 2 , 

4 2) = ^i 2) + (iP2-5i)ci + (2 g3 -3 g2 )c2, 

h< 2 } = H 2 2) - (5535i + 6<7i<?2 - 51P1 - §53£>2)ci - (-4^ - q\ + iq 2 q 3 + \q x y 2 + qip 3 )c 2 . 



They form two bi-Hamiltonian chains 



irodh^ = 

TT dh[ 1) = x[ 1} = mdh^ 

= TTld/l^ , 

with the corresponding Poisson operators ttq and 7Ti 



"0 = 



n 


1 












7Tl 



ir dh^ 

TTgdh^ 

■n^dh^ 



ni 



= 

= x[ 2) =n 1 dh { V 
= xP=n 1 dh[ 2) 

= 7Tld/4 2) 



X 



41) 



The separation curve for the extended system takes the form 

M(ciA + h^) + c 2 \ 2 + h {2) \ + hf = M 3 . 
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5 Summary 

We have considered the Stackel systems classified using their separation relations. The most general form 
of the separation relations considered in the present paper is 

m 

Y J Si{Vi,\)H {k \\ i )=il)i{\i,lJ,i), m<n, i = l,...,n, 
k=i 

where 

H ik) {^) = J2K k ~ J H ( 3 k \ n 1 + --- + n m = n 

and Sf,ipi are smooth functions of their arguments. Moreover, we have proved that all systems whose 
separation relations are of the above form admit (after the lift to an extended phase space) Gel'fand- 
Zakharevich bi-Hamiltonian representation. This confirms universality of the latter property for the 
Stackel systems. As a consequence, a geometric separability theory, based on the existence of GZ bi- 
Hamiltonian representation of a given system, is applicable for all Liouville intcgrable systems from the 
classes we considered. 

Quite obviously, the knowledge of quasi-bi-Hamiltonian representation is sufficient for separability 
of a given Liouville integrable system. Unfortunately, there is no systematic method available for the 
construction of such representation. On the other hand, there are some systematic methods for finding 
the bi-Hamiltonian representation. From this point of view the result presented in this paper is of interest, 
as it shows that the existence of bi-Hamiltonian representation is an inherent property of Stackel systems. 
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